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ABSTRACT : Inbreeding is the mating of relatives that produce progeny having more homozygous alleles than non-inbred animals. 
Inbreeding increases numbers of recessive alleles, which is often associated with decreased performance known as inbreeding 
depression. The magnitude of inbreeding depression depends on the level of inbreeding in the animal. Level of inbreeding is expressed 
by the inbreeding coefficient. One breeding goal in livestock is uniform productivity while maintaining acceptable inbreeding levels, 
especially keeping inbreeding less than 20%. However, in closed herds without the introduction of new genetic sources high levels of 
inbreeding over time are unavoidable. One method that increases selection response and minimizes inbreeding is selection of individuals 
by weighting estimated breeding values with average relationships among individuals. Optimum genetic contribution theory (OGC) uses 
relationships among individuals as weighting factors. The algorithm is as follows: i) Identify the individual having the best EBV; ii) 
Calculate average relationships (~) between selected and candidates; iii) Select the individual having the best EBV adjusted for average 

relationships using the weighting factor k, EBV* = EBV (l -kr. )• i v ) Repeat process until the number of individuals selected equals 

number required. The objective of this study was to compare simulated results based on OGC selection under different conditions over 
30 generations. Individuals (n = 1 10) were generated for the base population with pseudo random numbers of N~ (0, 3), ten were 
assumed male, and the remainder female. Each male was mated to ten females, and every female was assumed to have 5 progeny 
resulting in 500 individuals in the following generation. Results showed the OGC algorithm effectively controlled inbreeding and 
maintained consistent increases in selection response. Difference in breeding values between selection with OGC algorithm and by EBV 
only was 8%, however, rate of inbreeding was controlled by 47% after 20 generation. These results indicate that the OGC algorithm can 
be used effectively in long-term selection programs. (Key Words : Genetic Contribution Theory, Genetic Response, Breeding Value) 



INTRODUCTION 

Selection response has been improved continuously 
through advanced genetic selection programs recently 
although increase of inbreeding rate has been issued at the 
same time (Weigel and Lin, 2002). Inbreeding is the 
matings of individuals with ancestry in common and is 
based on the average relationship among all individuals in 
the population. Inbreeding increases homozygosity of 
alleles. Therefore, individuals with high inbreeding 
coefficients have a greater probability of homozygous 
alleles. Inbreeding has been used to develop uniform herds 
of animals with more homozygous genotypes, but increased 
inbreeding may cause side effects such as reproductive 
atrophy resulting in stillbirth or deformity. In pigs, 
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reduction of litter size due to increased inbreeding can be 
the result of either decreased conception rate as inbreeding 
level of sows increase, or loss of vitality of piglets as their 
inbreeding level increases. 

Falconer (1960) has reported on experiments from long- 
term, full-sib, mating of mice that approximately 40% of 
the reduction in litter size due to inbreeding was attributable 
to reduced fertility of the females and the remaining 60% to 
reduced survivability of the young. In cattle, it has been 
reported that a 1% increase in the inbreeding coefficient 
resulted in decrease of milk production by 37 kg and the 
production period by 13.1 days in Holstein cows (Smith et 
al., 1998). In pigs a 10% increase in inbreeding of young 
pigs has been reported to reduce litter size by 0.29 while a 
10% increase in inbreeding of the female results in the 
reduction of 0.40 head per litter (Johnson, 1990). In poultry, 
Sewalem et al. (1999) have reported that inbreeding 
depression was not so much for egg production traits, but an 
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increase of inbreeding was related to a delay in sexual 
maturity. It has been reported that 0.7% of inbreeding 
coefficients on average increased every year after 23 
generations without full- and half-sib matings, which 
resulted in 3 eggs decrease per year (Ameli et al., 1991). 

With this reason, many methods to maximize genetic 
progress minimizing the rate of inbreeding have been 
suggested and studied. Weigel (2001) has summarized the 
methods, which are restrictions on family size, creation of 
sublines, mating programs, modifications to BLUP, and 
optimum contribution theory. Optimal inbreeding theories 
suggest that breeding values would be maximized while 
minimizing inbreeding. It is said that methods based on 
optimum contribution theory are more promising among 
them (Konig et al., 2010). Although numerous cases were 
studied, there are few researches suggesting which 
constraints would be used properly in overall because they 
rely on specific breeding programs and the individual 
numbers of foundation population. One of the theories has 
been suggested by Quinton and Smith (1995), which is 
based on the consideration of relationship coefficients 
among selection candidates resulting in their re-ranking. It 
also depends on the size of the foundation herds. Their 
method is relatively intuitive than other methods in terms of 
the factors to be considered in the model such as the 
number of foundation group, heritability, and weighting 
factors. However, they just showed the results of 10 
generations, which may not be enough to explain the effect 
for building an inbred line. Also, it was not clear about the 
level of weighting factors in each gender group with 
different mating conditions. 

The purpose of this study is not to improve the existing 
method, nor to suggest a new algorithm. It is to evaluate the 
effect of OGC theory on the inbreeding through the 
computer simulation of the method suggested by Quinton 
and Smith (1995), and to find out how much it would 
control the rate of inbreeding in a closed herd starting 10 
sires and 100 dams with different combinations of 
weighting factors to breeding groups over 20 generations. 
Livestock producers will have some idea regarding to the 
optimal combinations of weighting factors in a particular 
mating condition from the results in this study. 

MATERIALS AND METHODS 

Increasing selection response with minimal inbreeding 
is based on controlling the ranking of individuals weighting 
selection based on EBVs with average relationships among 
individuals. Optimum genetic contribution theory (OGC) 
uses those relationships as weights. Selection based only on 
EBV from BLUP evaluations results in increased rates of 
inbreeding especially in long-term selection because the 
probability is greater that relatives with superior EBVs will 



be selected together. Optimum genetic contribution as 
suggested by Quinton and Smith (1995) controlled the 
increase of inbreeding in a long-term closed miniature pig 
selection program by considering EBV and the rate of 
inbreeding between candidates and other individuals in the 
population. The steps involved in the algorithm are as 
follows: 

i) Identify the individual having the best EBV; 

ii) Calculate average relationships (r.) between 

selected individual(s) and candidate individuals; 

iii) Select the individual having the best OGC score 
with EBV* = EBV j (l- ifcrj using the weighting factor k; 

iv) Repeat the steps until the number of individuals 
required are reached; 

v) Repeat process for the opposite sex. 

Quinton and Smith (1995) demonstrated that OGC 
could maintain genetic improvement comparable to 
selection on only EBVs in a simulation of 10 generations 
with k values from 0.5 to 2.0. In this study, 20 generations 
were simulated for 25 combinations of k for both males and 
females where k ranged from 0 to 2.0 and was increased by 
0.5. The simulated trait was litter size with a heritability of 
0.1. One hundred and ten pigs were generated for the base 
population with pseudo random numbers N~ (0, 3). Ten 
pigs were assumed as male, and the rest as female because 
most of populations show polygamic structures with high 
mating ratios (Sanchez et al., 2003). Each male was mated 
to ten females, and every female was assumed to have 5 
progeny resulting in 500 individuals for the following 
generations. Sex ratio was assumed to 1:1. Ten males and 
one hundred females were selected and mated using the 
OGC algorithm for the next generations. This same process 
of selection was repeated for 20 generations. This 
simulation was analyzed without considering the aspect of 
overlapping generations. Average EBVs and the rates of 
inbreeding after ten simulations were calculated each 
generation. 

RESULTS AND DISCUSSION 

These results show that the OGC algorithm simulated in 
this study can control inbreeding effectively while 
maintaining consistent selection response. Differences in 
breeding values between different matings with the OGC 
algorithm and with only EBVs were less than 8% (Tables 3 
and 4). However, the rate of inbreeding could be controlled 
as much as about 47% in this OGC simulation over 20 
generations (Tables 1 and 2). 

Table 3 and 4 are average breeding values and standard 
deviations over 20 generations of selection with different k 
values. The expected breeding values after 20 generations 
were about 8% higher with k = 0 than k = 2 for both sexes 
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Table 1. The percentage of increase in the rate of inbreeding with 
different k values for female when k for male = 0 



Generation 






k for female 






0.0 


0.5 


1.0 


1.5 


2.0 


1 


0.0 


0.0 


0.0 


0.0 


0.0 


2 


2.9 


2.6 


2.3 


1.9 


1.8 


5 


22.7 


19.9 


17.8 


16.7 


16.3 


10 


37.5 


26.5 


21.4 


20.6 


17.5 


15 


49.3 


28.2 


24.8 


21.8 


19.6 


20 


57.8 


31.9 


26.3 


23.4 


20.2 


Table 2. The percentage of increase in the rate of inbreeding with 


different k values for male when k for female = 


0 




Generation 






k for male 






0.0 


0.5 


1.0 


1.5 


2.0 


1 


0.0 


0.0 


0.0 


0.0 


0.0 


2 


3.1 


2.9 


2.7 


2.2 


1.9 


5 


22.7 


12.1 


9.9 


9.4 


7.6 


10 


37.5 


13.2 


10.5 


9.7 


8.0 


15 


49.3 


13.7 


10.9 


10.5 


10.3 


20 


57.8 


14.9 


12.5 


11.7 


10.6 



Table 3. Percent of increases in the rates of inbreeding compared 
to the case of k = 0.0 for male and female 



k k (female) 



(male) 




0.5 


1.0 


1.5 


2.0 


0.5 


5 


9.09 


7.82 


7.10 


5.78 




10 


10.38 


10.24 


9.77 


9.34 




15 


11.12 


10.83 


9.99 


9.87 




20 


13.30 


11.46 


11.41 


11.19 


1.0 


5 


8.80 


6.95 


6.93 


5.40 




10 


9.38 


8.18 


7.69 


7.26 




15 


9.57 


8.80 


8.33 


7.62 




20 


11.02 


9.42 


8.85 


8.81 


1.5 


5 


6.26 


5.26 


5.34 


4.34 




10 


7.22 


6.61 


6.21 


5.49 




15 


7.62 


7.39 


7.36 


5.95 




20 


9.42 


8.31 


7.88 


7.45 


2.0 


5 


4.89 


4.25 


4.12 


3.13 




10 


5.95 


5.49 


5.28 


4.37 




15 


7.10 


6.27 


5.60 


5.31 




20 


8.00 


7.12 


6.46 


5.98 



Table 4. Changes of average breeding values and standard deviations in male with different k values for female when k for male = 0 

k for female 

Generation 





0.0 


0.5 


1.0 


1.5 


2.0 


1 


3.00±0.04 


3.0110.14 


2.8510.03 


2.9910.17 


2.9910.03 


2 


3.65±0.15 


3.5910.10 


3.6710.07 


3.4110.14 


3.6710.05 


3 


4.1010.18 


4.0810.11 


4.1910.10 


3.9210.13 


4.1210.11 


4 


4.4910.09 


4.5810.19 


4.6310.03 


4.3410.05 


4.5010.12 


5 


4.9610.15 


4.9910.12 


5.0110.03 


4.7610.12 


4.8610.13 


6 


5.3810.12 


5.4110.20 


5.4210.01 


5.1210.11 


5.2510.09 


7 


5.7710.11 


5.8310.22 


5.8310.02 


5.5410.11 


5.6510.10 


8 


6.2110.10 


6.2610.24 


6.2410.01 


5.9410.13 


6.0510.11 


9 


6.6410.12 


6.7210.25 


6.6610.05 


6.3810.11 


6.4310.10 


10 


7.0710.11 


7.1310.23 


7.0810.02 


6.7610.13 


6.8210.13 


11 


7.4910.09 


7.5710.23 


7.4810.01 


7.1810.15 


7.2610.13 


12 


7.9210.13 


7.9810.14 


7.9510.01 


7.5610.15 


7.5810.09 


13 


8.3410.15 


8.4310.14 


8.3210.03 


7.9910.14 


7.9910.08 


14 


8.7610.13 


8.8210.16 


8.7210.06 


8.4110.15 


8.3410.08 


15 


9.2010.12 


9.2710.15 


9.1210.06 


8.8110.15 


8.6810.06 


16 


9.6310.11 


9.6610.16 


9.5610.04 


9.1610.08 


8.9310.08 


17 


10.0810.13 


10.0710.17 


9.9610.01 


9.5610.09 


9.3610.03 


18 


10.4810.10 


10.4710.19 


10.3610.03 


9.9210.09 


9.7010.15 


19 


10.9110.14 


10.8710.16 


10.8010.02 


10.29+0.05 


10.0410.13 


20 


11.3310.11 


11.3010.18 


11.2010.01 


10.6810.04 


10.4510.21 
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restricting k = 0 in one sex. However, overall breeding 
values were increased regardless of k values. Although 
inbreeding depression was not included in the simulated 
model, it should not be a problem to see the effect of the 
OGC algorithm in this study. 

Trends of inbreeding coefficients by generation based 
on male and female (Tables 1 and 2) selections show that 
the rates of inbreeding were controlled with different k 
values. As k values are higher, the rates of inbreeding are 
decreased. When both k values for male and female were 
zero, the rate of inbreeding of the 20 th generation was 
increased by 57.8% compared to the base generation; 
however, when k value was set to 2.0 in one sex, the rates of 
inbreeding were limited to 20.2% and 10.6% increase with 
k = 0.0 for male and for female, respectively, in the 20th 
generation. This is because higher k values emphasize the 
relationship between candidates, which results in changes of 
breeding values that would be lower than the original. 
Therefore, when k = 2.0 is compared to k = 1.0, more strict 
changes in the original ranks of candidates who have 
greater degrees of relationship occur than in the candidates 
having less degrees of relationship resulting in a less related 
group of candidates controlling the rate of inbreeding. In the 
mean time, the rates of inbreeding were just 9.4% increase 
for male selection and 8.3% increase for female selection 
with k = 2.0 in both sexes compared to the case of selection 
with k = 0.0 in both sexes (Table 3). 



From the results, the impact of k values in male on the 
level of inbreeding for the next generation was greater than 
that in the female because 10 of 250 were selected as male 
and 100 of 250 were selected as female assuming a sex ratio 
of 1:1 for the 500 simulated progeny. In other words, 
selection intensity was higher in male than in female 
population. Therefore, increasing the k value for the males 
had a larger effect on reducing inbreeding than increasing 
the k value of females. 

Methods of controlling inbreeding and the algorithm 
procedures have been well summarized by Weigel (2001). 
The research applying a cost factor to the average 
relationship in the selected group decreased the rate of 
inbreeding by 35 to 39% after 10 generations of selection 
with just 1 to 3% reduction of selection response (Wray and 
Goddard, 1994). This kind of results can be found in several 
researches. Meuwissen and Goddard (1998) showed that 
selection with adjusted EBV decreased a genetic response 
by 4% with young bulls for AI progeny test programs while 
reducing a mean relationship among selected individuals by 
22%. 

The results in this study indicate that the OGC 
algorithm can be used effectively for the long-term 
selection program of any closed population. Also, it shows 
the possibility that inbreeding could be easily controlled by 
weighting EBV of candidates in each sex group. Software 
to predict selection responses and rates of inbreeding from 



Table 5. Changes of average breeding values and standard deviations in female with different k values for male when k for female = 0 



Generation 






k for male 






0.0 


0.5 


1.0 


1.5 


2.0 


1 


0.51±0.04 


0.5910.04 


0.5310.02 


0.5010.15 


0.5610.11 


2 


1.84±0.04 


1.9310.01 


1.7910.15 


1.7710.05 


1.8110.02 


3 


2.59±0.12 


2.6110.07 


2.5510.14 


2.5110.02 


2.4810.09 


4 


3.16±0.15 


3.2010.05 


3.0010.16 


2.9310.02 


2.9410.12 


5 


3.62±0.14 


3.6910.01 


3.4610.15 


3.3710.01 


3.3510.06 


6 


4.12±0.15 


4.0910.01 


3.9610.17 


3.7710.01 


3.7310.06 


7 


4.57±0.13 


4.4610.03 


4.3410.19 


4.1610.02 


4.1410.02 


8 


4.99±0.10 


4.8310.02 


4.7510.18 


4.5810.01 


4.54+0.03 


9 


5.38±0.12 


5.2610.01 


5.1510.21 


5.0110.01 


4.9410.04 


10 


5.83±0.12 


5.6310.02 


5.5410.22 


5.4210.11 


5.42+0.03 


11 


6.27±0.09 


6.0510.01 


5.9210.22 


5.7810.05 


5.7610.03 


12 


6.70±0.10 


6.4310.03 


6.3510.20 


6.1910.10 


6.1310.01 


13 


7.11±0.13 


6.8110.05 


6.7710.22 


6.5910.07 


6.5310.05 


14 


7.55±0.14 


7.2110.06 


7.1010.22 


6.9510.06 


6.9810.04 


15 


7.9610.15 


7.6410.09 


7.4910.17 


7.3410.14 


7.3310.01 


16 


8.3710.16 


8.0310.06 


7.9010.14 


7.7310.15 


7.7310.01 


17 


8.8010.17 


8.4110.01 


8.3410.14 


8.1210.11 


8.1410.05 


18 


9.2510.15 


8.7810.02 


8.7210.15 


8.5010.11 


8.57+0.05 


19 


9.6610.11 


9.1810.01 


9.1410.13 


8.8810.12 


8.9610.02 
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breeding programs using this theory would be beneficial to 
the livestock producers as well as animal keepers. Rutten et 
al. (2002) have developed the software named as 
"SelAction" to show possible combinations for optimal 
mating structure. Although inbreeding depression of the 
trait was not included in the Mendelian random sampling of 
this simulation study, the results of this simulation would be 
valuable because the reduction of the rate of inbreeding was 
5 times more than the reduction of breeding values. 
Hinrichs et al. (2006) have developed an algorithm titled in 
OCSELECT that is more useful when the number of 
candidates is huge, reducing computing time. The 
relationship matrix is included in this algorithm as A = 
ZApZ v +D, where Ap is the relationship matrix of parents, D 
is a diagonal matrix having variances of Mendelian 
sampling, Z is genetic contributions from parents to 
progeny. 

The stimulation model in this study was not perfect 
because the inbreeding depression of the trait was not 
included in the Mendelian random sampling model. Many 
cases about inbreeding depression have been reported. Irish 
Holstein-Friesian dairy cows that were 12.5% inbred 
showed decreases of 61.8, 5.3, and 1.2 kg in milk, fat, and 
protein yields, respectively (Parland et al., 2007). This 
inbred group also was expected to have a 2% greater chance 
of dystocia, and a 1% greater chance of stillbirth as 
inbreeding may cause reproductive problems. Gulisija et al. 
(2007) have reported that the range of the rate of inbreeding 
was 0.6 to 34% in their study which was nonparametric 
analysis of the impact of inbreeding on production with 
59,778 Jersey cows. In their results, milk yields were 
affected by inbreeding depression when the rate of 
inbreeding was greater than 7%. Therefore, inbreeding 
depression in milk yields would be prevented with A>1.5 in 
both sexes from the results of this study (Table 3). 

It is obvious that OGC theory can reduce effectively the 
rate of inbreeding over generations with a minimal decrease 
of selection response. Therefore, the results of this study 
will help livestock producers focus on a plan to decide the 
starting numbers of foundation animals to construct a 
breeding herd. In addition, this kind of simulation studies 
will provide them their effectiveness and optimal conditions 
in each theory. Further simulation research would be needed 
to consider inbreeding depression in Mendelian sampling 
procedure as well as overlapping generations same as 
commercial selection programs in practice. 
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